Improved low Q 2 model for neutrino and electron nucleon cross sections 
in few GeV region 

Arie Bodek a , Inkyu Park a , and Un ki Yang b 

a Department of Physics and Astronomy, University of Rochester, Rochester, New York 14618, USA 
b Enrico Fermi Institute, University of Chicago, Chicago, USA 

We present an improved model on neutrino- and electron-nucleon scattering cross sections using effective 
leading order PDFs with a new scaling variable £ w . Non-perturbative QCD effects at low Q 2 are separately 
treated for u and d valence quarks, and sea quarks. The improved model uses all inelastic charged lepton Fi data 
(SLAC/BCDMS/NMC/HERA), and photoproduction data on proton and deuterium. In this way, we obtain an 
improved model which describes all inelastic scattering charged lepton data, JLAB resonance data, and neutrino 
data. This improved model is expected to be better for neutrino oscillation experiments at few GeV region. 
(Presented by Un Ki Yang at Nulnt04, Mar. 2004, Laboratori Nazionali del Gran Sasso -INFN - Assergi, Italy) 



Recent discovery of neutrino oscillation from 
atmospheric neutrino experiment U have indi- 
cated that it is important to understand neutrino 
cross section in the few GeV region. But we do 
not have precise neutrino data to understand this 
region very well. Thus, good modeling of neutrino 
cross sections at low energies becomes very cru- 
cial for the precise next generation neutrino os- 
cillations experiments. In few GeV region, there 
are three types of neutrino interactions, quasi- 
elastic, resonance, and inelastic scattering. It 
is very challenging to disentangle each contribu- 
tion separately, especially, resonance scattering vs 
deep inelastic scattering (DIS) contribution, due 
to large non-perturbative QCD corrections to the 
DIS contribution. 

Our approach is that we describe these pro- 
cesses in terms of quark-parton model using 
the precise charged-lepton scattering data. It 
is then simpler to convert charged-lepton scat- 
tering cross section into neutrino cross sec- 
tion. In a previous report 0, we showed 
that our effective leading order model us- 
ing an improved scaling variable £ w describes 
all deep inelastic scattering charged lepton- 
nucleon scattering data including resonance data 
(SLAC/BCDMS/NMC/HERA/JLab) @EJ from 
very high Q 2 to very low Q 2 (down to photo- 



production region), as well as CCFR neutrino 
data citeyangthesis,rccfr. 

Our proposed scaling variable, £ w is derived on 
following basis. Using energy momentum conser- 
vation, the factional momentum, £ carried by a 
quark in a proton target of mass M can be ob- 
tained as follows; 

f = *Z41 (i) 

Q 2 (l + Vl + (2Mx) 2 /Q 2 )' 

where 

2Q' 2 = [Q 2 + M f 2 - M 2 } 

+ ^{Q 2 + M f 2 - M % 2 ) 2 + AQ 2 {KU 2 + P 2 ) 

Here Mi is the initial quark mass with average 
initial transverse momentum Pt and Mf is the 
mass of the quark in the final state. This expres- 
sion for £ was previously derived [B] for the case 
of Pt — 0. Assuming Mj = we use following 
variable in our model; 

_ 2x(Q 2 +M f 2 + B) 

Q 2 [l + ^l + {2Mx) 2 /Q 2 ]+2Ax' 

here, Mf = except charm-production case 
(Mf =1.5 GeV) in neutrino scattering. The pa- 
rameter A accounts for the higher order (dynamic 
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higher twist) QCD terms in the form of an en- 
hanced target mass term (the effects of the proton 
target mass are already taken into account using 
the exact form in the denominator of t; w ) . The 
parameter B accounts for the initial state quark 
transverse momentum and final state quark ef- 
fective AM/ 2 (originating from multi-gluon emis- 
sion by quarks). This parameter allows that we 
could describe the photoproduction limit (all the 
way down to Q 2 =0). 

A brief summary of our effective leading order 
(LO) model is given as follows; 

1. The GRV98 LO PDFs 9] are used to de- 
scribe the F 2 data at high Q 2 region: The 
minimum Q 2 value of this PDFs is 0.8 
GeV 2 . 

2. The scaling variable x is replaced with the 
improved scaling variable £ w in Eq. [21 

3. Like earlier non-QCD based fits ^U] to low 
energy data, we multiply all PDFs by K 
factors, depending on the types of PDFs; 

K sea (Q 2 ) = Q2 Q +Cs 

K va lence{Q 2 ) = [1 ~ G D (Q)) 

where G D = 1/(1 + Q 2 /0.71) 2 is the proton 
elastic form factor. At low Q 2 , [1 — G 2 D (Q 2 )] 
is approximately Q 2 / (Q 2 + 0.178), which is 
very close to our earlier fit result [3] . These 
modifications were done in order to describe 
low Q 2 data in the photoproduction limit, 
where F 2 is related to the photoproduction 
cross section according to 

47r 2 a E M „ 0.112m& . . 

a ^P) = Q2 F 2 = Q2 F 2 (4) 

4. We freeze the evolution of the GRV98 PDFs 
at a value of Q 2 = 0.80. Below this Q 2 , F 2 
is given by; 

F 2 {x, Q 2 < 0.8) = K{Q 2 ) x F 2 (£, Q 2 = 8) (5) 



5. Finally, we fit to the effective GRV98 
LO PDFs (£ w ) using the DIS data 
(SLAC/BCDMS/NMC/H1). We ob- 
tained excellent fitting results; ^4=0.419, 
5=0.223, and 0,1=0.544, 0,2=0.431, and 
C sea =0.380 and X 2 /DOF = 1235/1200. 
Because of the K factors to the PDFs, we 
find that the GRV98 PDFs need to be in- 
creased by 7V=1.011. 

In the previous analysis, the structure functions 
data were corrected for the BCDMS systematic 
error shift and for the relative normalizations be- 
tween the SLAC, BCDMS, NMC and HI data. 
The deuterium data were corrected for nuclear 
binding effects |11I12| . We also included charm 
production contribution using the photon-gluon 
fusion model in order to fit the very high v HERA 
data. This was not necessary for any of the low 
energy comparisons but was only required to de- 
scribe the highest v HERA F3 and photoproduc- 
tion data, since the GRV98 PDFs did not include 
the charm sea for Q 2 > 0.8 GeV 2 . 

Performance of this effective LO model was 
very good in describing various DIS data from 
high Q 2 to low Q 2 region. Predictions for photo- 
production data on proton and deuteron (Q 2 = 
limit) showed good agreement. Furthermore, this 
model showed a reasonable description of the av- 
erage value of F 2 for SLAC and Jefferson reso- 
nance data [5] (Q 2 down to 0.07). 

In this report, we improve our effective LO 
model by treating low Q 2 corrections (K fac- 
tor in Eq. EJ) for u, d valence and sea quarks 
separately. Since the predictions of our previ- 
ous model showed good agreements with pho- 
toproduction data, we now include these data 
to get a better constraint on our model, in- 
stead of predicting. The fitted results using the 
improved model are given as follow: ^4=0. 538, 
5=0.305, CW=0.202, C olu =0.291, (7^=0.255, 
a 2u =0.189, C s ld=0.621, C s i u =0.363, and 
fPDF=lMh, and x 2 /DOF =1874/1574. 

Comparison of all DIS F 2 data and our fits are 
shown in Figure^and Our effective LO model 
describes the DIS F 2 data at low x and high x 
regions as well. This model also fits the photo- 
production data on proton and deuteron targets 
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very well, as shown in Figure [3] 

Comparisons of predictions using our model to 
the resonance data [5] which are not included in 
the fit is shown in Figure 0] Based on duality 
argument [15) . our model provides a reasonable 
description of the average value of F 2 data in the 
resonance region (down to Q 2 = 0.07). 

In addition to vector structure function, we 
have an axial vector structure function in neu- 
trino scattering. At the Q 2 — limit, the vector 
structure function should go to zero, while the 
axial-vector part has a finite contribution. At 
high Q 2 , these two structure functions are ex- 
pected to be same. Thus, it is important to 
understand the axial vector contribution at low 
Q 2 using the low energy neutrino data. As a 
preliminary step, we compare the CCFR j(i!7j 
and CDHSW energy neutrino data with our 
model, assuming that the vector contribution 
is same as the axial vector contribution. Fig- 
ure |S] show comparisons of our model to the high 
energy CCFR and CDHSW neutrino charged- 
current cross section data. In this comparison we 
correct for nuclear effects in iron. The structure 
function 2xF\ is obtained by using the R WO rid 
fit from reference 0]. Our predictions show very 
good agreement with these neutrino data on iron. 

In order to have a full description of all charged 
current processes, the contribution from quasi- 
elastic scattering must be added separately at 
x = 1. The best prescription is to use our 
model in the region above the first resonance 
(above W—1.35 GeV) and add the contributions 
from quasi-elastic and first resonance (W=1.23 
GeV) separately. This is because the W=M 
and 1^=1.23 GeV regions are dominated by one 
and two isospin states, and the amplitudes for 
neutrino versus electron scattering are related 
via Clebsch-Gordon rules ^2 instead of quark 
charges (also the V and A couplings are not equal 
at low W and Q 2 )- In the region of higher mass 
resonances (e.g. W—1.7 GeV) there is a signifi- 
cant contribution from the deep-inelastic contin- 
uum which is not well modeled by the existing 
fits ^3] to neutrino resonance data (and using 
our modified PDFs should be better). For nuclear 
targets, nuclear corrections [3j must also be ap- 
plied. Recent results from Jlab indicate that the 
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Figure 1. Our fit results to the effective LO 
model is compared to the F 2 data (SLAC, 
BCDMS, NMC) at high x:[top] F 2 proton, [bot] 
F 2 deuteron. 
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Figure 2. Our fit results to the effective LO model 
are compared to the F% data at low x: [top] NMC 
F 2 data, [bot] HI F 2 data. 
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Figure 3. Our fit results to the effective LO model 
are compared to the photoproduction data (Q 2 = 
limit) [top] proton, [bot] deuteron. The dotted 
line do not include charm contribution from gluon 
fusion process 
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Fe/D ratio in the resonance region is the same as 
the Fe/D ratio from DIS data for the same value 
of £ (or £ w ). 

We are currently working on constraining the 
low Q 2 axial vector contribution using low energy 
CDHSW and CHORUS [TB] data. The form of 
the fits we plan to use is motivated by the Adler 
sum rule |17) for the axial vector contribution as 
follows: 

Q 2 
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2v — ax 



Q + Civ-c 



(6) 



where F A (Q 2 ) = -1, 267/(1 + Q 2 /1.00) 2 . 

In addition, we plan to implement a modulat- 
ing function |18I19| A(W,Q 2 ) to improve model- 
ing in the resonance region by including (instead 
of predicting) the resonance data. Although a de- 
scription of the average cross section in the res- 
onance region is sufficient for most neutrino ex- 
periments, because of the effects of experimental 
resolution and Fermi motion for nuclear targets. 
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Figure 4. Comparisons of the resonance data and 
the predictions of our effective LO model: [top] 
proton, [bot] deuteron. 



Figure 5. [top] Comparisons of the CCFR neu- 
trino data and the predictions (solid line) of our 
effective LO model at 55 GcV, the dotted line us- 
ing the default GRV98 PDFs: [bot] Ratio of the 
CDHSW neutrino data and the predictions of our 
effective LO model at 23 GeV. 



